Efficient Quantum Computing of Complex Dynamics 
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We propose a quantum algorithm which uses the number of qubits in an optimal way and 
efficiently simulates a physical model with rich and complex dynamics described by the quantum 
sawtooth map. The numerical study of the effect of static imperfections in the quantum computer 
hardware shows that the main elements of the phase space structures are accurately reproduced 
up to a time scale which is polynomial in the number of qubits. The errors generated by these 
imperfections are more dangerous than the errors of random noise in gate operations. 
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When applied to computation, quantum mechanics 
opens completely new perspectives: a quantum com- 
puter, if constructed, could perform certain computations 
faster than classical computers, exploiting quantum me- 
chanical features such as entanglement or superposition 
m. Shor discovered a quantum algorithm 0] which fac- 
torizes large integers exponentially faster than any known 
classical algorithm. It was also shown by Grover H that 
the search of an item in an unstructured database can 
be done with a square root speed up over any classi- 
cal algorithm. However, at present only a few quantum 
algorithms are known which simulate physical systems 
with exponential efficiency. Such systems include certain 
many-body problems Q, spin lattices 1^, and models of 
quantum chaos 



In this Letter, we present a quantum algorithm which 
computes the time evolution of the quantum sawtooth 
map, exponentially faster than any classical computa- 
tion. This model has rich and complex dynamics and 
finds various applications, e.g. for dynamical localiza- 
tion in billiards |q-[lC|]. Our algorithm, based on the 
Quantum Fourier Transform (QFT) [^, simulates the 
dynamics of a system with N levels in 0((log2 N)'^) op- 
erations per map iteration, while a classical computer, 
which performs Fast Fourier Transforms (FFT), requires 
0(A^log2 A^) operations. A further striking advantage of 
the algorithm is the optimum utilization of qubits: one 
needs only Uq = log2 N qubits (without any extra work 
space) . We demonstrate that complex phase space struc- 
tures can be simulated with less that 10 qubits, while 
about 40 qubits would allow one to make computations 
inaccessible to present-day supercomputers. This is par- 
ticularly important, since experiments with few qubits 
are being performed at present [ p2||l3| , in particular the 
QFT was implemented on a three qubit nuclear mag- 
netic resonance quantum computer JMJ. For this reason 
the investigation of this interesting physical system will 
be accessible to first quantum computers, operating with 



few qubits and for which large-scale computations like 
integer factoring are not possible. 

The classical sawtooth map is given by 



n = n + k{9 — tt) , 



Tn, 



(1) 



where (n, 9) are conjugated action-angle variables (0 < 
6 < 27r) , and the bars denote the variables after one map 
iteration. Introducing the rescaled variable p = Tn, one 
can see that the classical dynamics depends only on the 
single parameter K = kT, so that the motion is stable 
for —4 < K < and completely chaotic for K < —4 and 
K > 0. For such a discontinuous map the conditions of 
the Kolmogorov-Arnold-Moser (KAM) theorem are not 
satisfied and therefore, for any K ^ 0, the motion is not 
bounded by KAM tori. The map (||) can be studied on 
the cylinder (p e (— oo, -|-oo)), which can also be closed to 
form a torus of length 2ttL, where L is an integer. For any 
K > 0, one has normal diffusion: < (Ap)^ >w D{K)t, 
where t is the discrete time measured in units of map 
iterations and the average < • • • > is performed over an 
ensemble of particles with initial momentum pq and ran- 
dom phases < 9 < 27r. It is possible to distinguish two 
different dynamical regimes Il3| : ior K > 1, the diffusion 
coefficient is well approximated by the random phase ap- 
proximation, D{K) w {Tr'^/3)K^, while for < X < 1 
diffusion is slowed down, D{K) w i.iK^^^, due to the 
sticking of trajectories close to broken tori (cantori). For 
—4 < K < {) the motion is stable, the phase space has 
a complex structure of elliptic islands down to smaller 
and smaller scales, and we observed anomalous diffusion, 
< (Ap)2 >(x i", (for example, a = 0.57 when i^ = -0.1). 
The quantum evolution on one map iteration is de- 
scribed by a unitary operator U acting on the wave func- 
tion Tp: 

" ^^'/V, (2) 

1). The classical limit 
corresponds to fc — > oo, T ^ 0, and K = kT = const. In 
this quantum model one can observe important physical 
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where n = —id/ 89 (we set h 



phenomena like dynamical localization ||,^. Indeed, due 
to quantum interference effects, the chaotic diffusion in 
momentum is suppressed, in a way similar to Anderson 
localization in disordered solids. Also in the vicinity of a 
broken KAM torus, cantori localization takes place, since 
a cantorus starts to act as a perfect barrier to quantum 
wave packet evolution, if the flux through it becomes less 
than^i |l|,|l^,|-|lO|. 

The most efficient way to simulate the quantum dy- 
namics (0) on a classical computer is based on for- 
ward/backward FFT between 6 and n representations. 
This is advantageous because the evolution operator JJ 
is the product of two unitary operators, C/fc = exp(ifc(0 — 
7r)V2) (kick) and Ut = exp(-JTnV2) (free rotation), 
which are diagonal in the 9 and n representation, respec- 
tively. Therefore, for a system with N levels, the one 
map iteration ^ requires two FFT and two diagonal 
multiplications and can be performed in 0{N\og2{N)) 
operations. The dynamics (g) can be simulated expo- 
nentially faster on a quantum computer with nq = log2 N 
qubits by means of the following quantum algorithm: (i) 
the wave function |V') = X]n=o ^"l") (given in the n 
representation) is multiphed by Ut, so that C/tIV") = 
^jja„exp(— iTn^/2)|n). This step can be done in n^ 
controUed-phase shift gates, as explained in 0; (ii) one 
can get the wave function in the 9 representation via 
the QFT |ll| , which requires Uq singlc-qubit (Hadamard) 
gates and nq{nq — l)/2 two-qubit gates (controlled phase- 
shifts) ; (iii) the action of Uk is diagonal in the angle rep- 
resentation and can be simulated in a way similar to (i) in 
ni two-qubit gates (we note that this is possible thanks 
to the particular form of Uk for the sawtooth map) ; (iv) 
we go back to the momentum basis performing backward 
QFT in nq{nq + I)/2 gates. Therefore the whole algo- 
rithm requires Ug = 3nl + Uq quantum gates per map 
iteration. 

We model the quantum computer hardware as a two- 
dimensional lattice of qubits (spin halves) with static 
fluctuations/imperfections in the individual qubit ener- 
gies and residual short-range inter-qubit couplings. The 
model is described by the following many-body Hamilto- 
nian (see also @): 



Hs 






{Ao + S,)a^ 



2^ Jjjcrfo-J 



(3) 



where the Ci are the Pauli matrices for the qubit i, and 
Aq is the average level spacing for one qubit. The sec- 
ond sum in (0) runs over nearest- neighbor qubit pairs, 
and Si, Jij are randomly and uniformly distributed in 
the intervals [—5/2,5/2] and [—J,J], respectively. We 
study numerically the many-body dynamics of the quan- 
tum computer (ra) running the quantum algorithm de- 
scribed above. The algorithm is realized by a sequence of 
instantaneous and perfect one- and two-qubit gates, sepa- 
rated by a time interval Tg , during which the Hamiltonian 



(^ gives unwanted phase rotations and qubit couplings. 
A special one-qubit rotation is applied to eliminate the 
average phase accumulation given by Aq. We consider 
the case of short range inter-qubit interaction, in which 
two-qubit gates can be performed only between nearest 
neighbor pairs. As a result, for a square lattice of Uq 

1/2 

qubits, we make 0{nq ) swap gates to obtain a generic 
two-qubit gate. Therefore a single iteration of the map 
(g) requires Ug = 0{nq' ) gates. 

In this Letter, we study the sawtooth map in the 
anomalous diffusion regime, with K = —0.1, — tt < p < tt 
(torus geometry). The classical limit is obtained by in- 
creasing the number of qubits Uq, with T = 2Tr/N (k — 
K/T, -N/2 < n < N/2). We consider as initial state 
at time i = a momentum eigenstate, |V'(0)) — l"o)i 
with no = [0.38A^]. Such a state can be prepared in 
0{nq) one-qubit rotations starting from the ground state 
|0, . . . , 0). The dynamics of the sawtooth map reveals 
the complexity of the phase space structure, as shown 
by the Husimi functions of Fig.|l|, taken after 1000 map 
iterations. We note that riq = 6 qubits are sufficient to 
observe the quantum localization of the anomalous diffu- 
sive propagation through hierarchical integrable islands. 
At Uq = 9 one can see the appearance of integrable is- 
lands, and at nq = 16 the quantum Husimi function ex- 
plores the complex hierarchical structure of the classical 
phase space. The effect of static imperfections (0) for the 
operability of the quantum computer is shown in Fig.n^ 
(right column). The data are shown for J = (we ob- 
served similar structures for J = 5). The main features 
of the wave packet dynamics remain evident even in the 
presence of significant imperfections, characterized by the 
dimensionless strength e = 5Tg. The main manifestation 
of imperfections is the injection of quantum probabil- 
ity inside integrable islands. This creates characteristic 
concentric ellipses, which follow classical periodic orbits 
moving inside integrable islands. These structures be- 
come more and more pronounced with the increase of 
Uq. Thus quantum errors strongly affect the quantum 
tunneling inside integrable islands, which in a pure sys- 
tem drops exponentially (ex exp(— CiV), C =const). It 
is interesting to stress that the effect of quantum errors 
is qualitatively different from the classical round-off er- 
rors, which produce only slow diffusive spreading inside 
integrable islands (see Fig.|l|). This difference is related 
to the fact that spin flips in quantum computation can 
make direct transfer of probability on a large distance in 
phase space. 

A more quantitative study of the effect of static imper- 
fections is obtained from the fidelity of quantum compu- 
tation, defined by /(i) = |('0e(OlV'o(i))P, where ■(/'c(i) is 
the actual quantum wave function in the presence of im- 
perfections and ipo{t) is the quantum state for a perfect 
computation. The data of Fig.||show that f{t) drops 




FIG. 1. Husimi function for the sawtooth map in ac- 
tion angle variables {p,6), with — tt < p < n (vertical 
axis) and < ^ < 27r (horizontal axis), for K = —0.1, 
T = 27r/2"«, no = po/T = [0.38 x 2"«], averaged in the 
interval 950 < t < 1000. From top to bottom: n, = 6, 9, 16 
and classical density plot, obtained from an ensemble of 10 
trajectories, with initial momentum po = 0.38 x 27r and ran- 
dom angles. Left and right columns show the case without 
and with imperfections: in the quantum case the imperfec- 
tion strength e = Srg scales oc n~^ , where e = 2 x 10"'^ 
(rig = 6), e = 6 x 10"* (n, = 9), e = 10""' (n, = 16), at 
J = 0; in the classical case round-off errors are of amplitude 
10~^. We choose the ratio of the action-angle uncertainties 
s = Ap/Ae = 1 {ApAe = r/2). The color is proportional to 
the density: blue for zero and red for maximal density. 



with t according to the law f{t) w exp(— Ai^). This 
can be understood as follows: the fluctuations in the 
individual qubit energies, coupled to the gate opera- 
tions in time, give effective Rabi oscillations. For this 
reason f{t) ^ \Yl^liCOs{SiTgngt)\'^ and at short times 
f{t) ^ exp(— nq(engi)^). This applies also for the case 
with quantum chaos in the Hamiltonian (0) (e.g. J = S), 
which at short time scale can be considered from the 
viewpoint of Rabi oscillations. Indeed, here we con- 
sider small e and Uq values, for which quantum chaos 
develops on a chaotic time scale t^ ~ l/J^/n^ » Tg 
p8[ . For example, for Uq — 9, e — 10"*, and J — 5 
one has t^ « 3 x lO^Tg, while one map iteration takes 



time Ti 



= 413Tg w r^^/lO. The fidelity time 



scale tf can be obtained from the condition /(i/) — c, 
for example c = 0.9. The above argument gives tf ex 
e^^Tiq^Uq (x e~^n~^. This estimate is in a good agree- 
ment with the numerical data shown in Fig.pl. These data 
also show that tf is not strongly affected by the fact that 
the internal quantum computer hardware (j^ is in the 
integrable (J = 0) or quantum chaos (J = S) regime. 
The most important point is that the dependence of the 
fidelity time tf on Uq is only polynomial. 
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FIG. 2. Fidelity as a function of t, for Uq — 9 qubits, J — 5 
(solid lines) and J — (dashed lines). From top to bottom: 



e= 10"^3x 10"^ 



10"*, 3 x 10"*, 10" 



It is interesting to compare the effect of static imperfec- 
tions generated by the Hamiltonian (0) in the regime of 
perfect gates with the case when the static imperfections 
are absent but the computer operates with noisy gates. 
To model the noisy gates we put J — and Si fiuctuating 
randomly and independently from one gate to another in 
the interval [—6/2,6/2]. In this case the fidelity drops 
in a qualitatively different way f{t) « exp{—Bt) (data 
not shown). This corresponds to the Fermi golden rule 
regime, where the probability at a given state decays ex- 
ponentially with time. At each gate operation a proba- 



bility of order e^ is transferred to other states, so that 
the fidehty time scale is given by tf ^ l/{e^ng). This es- 
timate is in good agreement with numerical data shown 
in Fig.^ We found the same behavior ior J — 5 — and 
noisy gates with unitary rotations on a random angle e. 
The dependence on e is qualitatively different comparing 
to the static imperfection case (||) discussed above. We 
stress that static imperfections give shorter times tf and 
therefore are more dangerous for quantum computation. 
We note that the quadratic and linear regimes of fidelity 
decrease with time were respectively discussed in Refs. 

PI,po|. 



point. This can also provide important new informations 
about quantum states in systems with hierarchical phase 
space structures. We believe that in the near future our 
results can be experimentally observed in quantum com- 
puters operating with a small number of qubits. 
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FIG. 3. Fidelity time scale tf as a function of e, for n, = 9, 
in the case of static imperfections {J — 5 (circles) and J — 
(squares)) and noisy gates (diamonds, see text). The straight 
lines have slopes —1 and —2. The inset shows the depen- 
dence of tf on the number of qubits, for e — W~*, J = 0; 
the power-law fit (straight line) gives tf oc n^^'®. Error bars 
give the statistical errors obtained from 10 < fe < 500 random 
configurations of static imperfections. 



In summary, we have shown that a realistic quantum 
computer can simulate with exponential efficiency com- 
plex dynamics with rich phase space structures. Our 
studies demonstrate that the main structures inside the 
localization domain are rather stable in the presence of 
static imperfections. We show that the errors from static 
imperfections are stronger than the errors of noisy gates. 
Of course it is not possible to extract all exponentially 
large information hidden in the wave function with 2"' 
states. However, it is possible to have access to coarse 
grained information. For example from a polynomial 
number of measurements one can obtain the probabil- 
ity distribution over momentum (or angle) states. This 
allows one to study the anomalous diffusion in the deep 
semiclassical regime. Such an information is not accessi- 
ble for classical computers which cannot simulate more 
than 2^" quantum states. Recently an efficient algorithm 
was proposed in Ref. pit], which allows one to measure 
the value of the Wigner function at a chosen phase space 
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